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The thermodynamics of the antiferromagnetic ordering transition in NiO and the photoelectron 
spectra in the antiferromagnetic phase are studied by the Variational Cluster Approximation. Using 
realistic Racah parameters to describe the Coulomb interaction in the Ni 3d shell and a Slater-Koster 
parameter (pda) which is slightly (10%) increased over the band structure estimate the calculated 
Neel temperature is 481 Kelvin (experimental value: 523 Kelvin). The magnetic susceptibility above 
Tn has Curie-Weiss form. A significant contribution to the stabilization of the antiferromagnetic 
phase comes from electron hopping between oxygen which would be missed in theories that consider 
superexchange along a single bond only. The single particle spectral function in the ordered phase 
is in good agreement with experiment, in particular a number of dispersionless bands which are not 
reproduced by most calculations are obtained correctly. These flat bands are shown to be direct 
experimental evidence for a dispersionless electronic self-energy with several poles in the energy 
range of the valence band which originate from the multiplets of the Ni 3+ ion. Small but possibly 
experimentally detectable changes of the photoelectron spectra with temperature are discussed, in 
particular a widening of the insulating gap in the paramagnetic phase by approximately 10% is 
predicted. 

PACS numbers: 74.20.Mn,74.25.Dw 


I. INTRODUCTION 

Nickel Oxide has received attention over several 
decades because it is the prototype of a correlated in¬ 
sulator. As early as 1937 this material was cited as a 
counterexample to the Bloch theory of solidsjlj: assum¬ 
ing strong ionicity, oxygen will be O 2- leaving Nickel to 
be Ni 2+ or [Ar]3d 8 . This means that 4 electrons/spin 
direction have to be distributed over five Ni 3d bands 
which could result in an insulating ground state only if 
one of the five bands were splitt off from the others over 
the entire Brillouin zone. This is impossible, however, 
because at T the five Ni 3d bands converge into one 3- 
fold degenerate (t- 2 g ) and one 2-folcl degenerate (e g ) level. 
Band structure calculations for the paramagnetic phase 
of NiO 3 confirm this, showing three O 2p derived bands 
well below the Fermi energy and a group of five Ni 3d 
bands which are intersected by the Fermi energy. 

This simple picture is modified in that NiO undergoes 
an antiferromagnetic ordering transition at Tjv = 523 
Kelvin. Thereby Ni-ions in planes perpendicular to 
(1,1,1) allign their magnetic moments parallel to each 
other, with the ordered moment in successive planes be¬ 
ing antiparallel (type II antiferromagnetism). This struc¬ 
ture is consistent with the Goodenough-Kanamori rules 
because all 180° Ni-O-Ni bonds are antiferromagnetic in 
this way. Band structure calculations within the frame¬ 
work of Density Functional Theory (DFT) for the antifer¬ 
romagnetic phase reproduce the insulating ground state, 
but the band gap is only G ~ 0.3 eU[s] whereas the ex¬ 
perimental value is G ~ 4.3 eV"Q. It should be noted, 
however, that DFT does indeed give a rather accurate 
estimate of G = 4.1 eV for the single-particle gap in 
NiO Q if one does not use the band structure of Kohn- 


Sham eigenvalues - which have no true physical signifi¬ 
cance anyway - but calculates the ground state energy 
Eq of finite clusters as a function of electron number N 
and uses G = E 0 (N + 1) + E 0 (N - 1) - 2 E 0 (N). The 
crucial point is, however, that NiO remains an insulator 
even above X)v so that the insulating nature of NiO can¬ 
not be explained by antiferromagnetic ordering. 

There is general agreement by now that the true origin 
of the insulating nature of NiO is the strong Coulomb 
interaction between electrons in the Ni 3d shell so that 
an adequate description requires a more accurate treate- 
ment of the electron-electron interaction. Accordingly, 
a wide variety of methods for treating interacting elec¬ 
trons have been applied to NiO over the years. Amongst 
others there are calculations using the self-interaction 
corrected density functional theory|gj, 0 ], the LDA+U 
formalism [1,0] and the GW approximation (Tol ITU ] which 
more recently, was also combined with the LDA+U 
formalism[H|. NiO was also treated in the frame¬ 
work of the three-body scattering formalism [+1 14) and 
dynamical mean-field theory (DMFT), both for the 
paramagnetic fl5l - fl7| and antiferromagnetic fl8j phase. 
Moreover it was shown recently that within DFT the 
agreement between calculated and measured band gap 
for antiferromagnetic NiO is improved considerably if a 
more accurate density functional is used ( 13 - 
A quite different - but very successful - approach was 
initiated by Fujimori and Minami[20j|. These authors 
showed that good agreement between theory and experi¬ 
ment could be obtained for angle-integrated valence band 
photoemission spectra if one focused on local physics by 
considering an octahedron-shaped NiOg°~ cluster com¬ 
prising of a single Ni-ion and its six nearest neighbor 
oxygen ions. The eigenstates and eigenenergies of such a 
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finite cluster can be calculated exactly by the configura¬ 
tion interaction (Cl) (or exact diagonalization) method 
and the single-particle spectral function be obtained from 
its Lehmann representation. The Cl method was subse¬ 
quently applied to the calculation of the angle-integrated 
valence band photoemission_spectra of a number of tran¬ 
sition metal compounds [21 h 29| | and was extended to sim¬ 
ulate X-ray absorption spectra [dOl - ldSL In all cases the 
agreement with experiment is excellent and in the case 
of X-ray absorption spectroscopy the comparison of sim¬ 
ulated and measured spectra by now is in fact becoming 
a routine tool for determining the valence and spin state 
of transition metal ions in solids [39l l40| . 

The reason for the success of the cluster method is that it 
takes into account the full Coulomb interaction between 
electrons in the transition metal 3d shell in the frame¬ 
work of atomic multiplet theory [4ll - [43l |. Introducing the 
compound index v = (n, l, to, a) (where n = 3 and l = 2 
for a 3d-shell) the Coulomb interaction between electrons 
in atomic shells can be written as [4l| - |43l | 

H 1 = \ V iyi’ v 3’ u k,n) ct.cl.C vl C„ k , 

i,j,k,l 

V(f 1 ,U2,l'3,l'i) = 6 cri,<T4 ^(72,CT3 ^mi+?722,m3-{-m4 

y: F k c fe (liTOi; Z4TO4) c k (l 3 m 3 \l 2 ni2)■ 
k 

(1) 

Here the c k (lm3',lm 2 ) denote Gaunt coefficients the F k 
Slater integrals and for a d-shell the multipole index 
k £ {0,2,4}. The Hamiltonian (JTJ) was derived orginally 
to explain the line spectra of atoms and ions in vacuum, 
see Ref. (4l| for an extensive list of examples. There 
are various simplified expressions in the literature [H| 
where the Hamiltonian © is approximated in terms 
of Hubbard-!/ and Hund’s rule J and sometimes addi¬ 
tional parameters (see Ref. 01 for the relation between 
Hubbard-!/ and Hund’s rule J and the Slater integrals) 
but the Hamiltonian © is the only one that can actually 
be derived from first principles and gives correct results 
for free ions (an instructive comparison of the eigenvalue 
spectra of the original Hamiltonian (HD and various sim¬ 
plified versions was given by Haverkort 0]). Hi contains 
diagonal terms such as 

( V(i/ 1 ,v 2 ,i'i,v2) - V(u 1 ,i/2,v 2 ,i'i) ) 1 Tlis 2 ( 2 ) 

but also off-diagonal terms where all four Vi in © are 
pairwise different. The off-diagonal matrix elements are 
frequently discarded in DMFT calculations flEl [Til be¬ 
cause they exchange electrons and thus exacerbate the 
minus-sign problem in quantum Monte-Carlo calcula¬ 
tions. On the other hand the matrix elements of these 
terms are of the same order of magnitude (namely pro¬ 
portional to the Slater integrals F 2 and F 4 ) as the dif¬ 
ferences between the various diagonal matrix elements in 


© so that there is no justification for discarding them 
but keeping different diagonal matrix elements. 

As will be seen below the Variational Cluster Approxi¬ 
mation (VCA) proposed by Potthoff [47l - l49| allows to ex¬ 
tend the scope of the Cl method of Fujimori and Minami 
once more, to the calculation of thermodynamical quan¬ 
tities and band structures for strongly correlated electron 
systems. Since the VCA is based on exact diagonaliza¬ 
tion and therefore free from the minus-sign problem the 
full Coulomb Hamiltonian © including the off-diagonal 
matrix elements can be included. As could have been 
expected on the basis of the considerable success of the 
Cl method in reproducing experimental spectra 0-0 
the VCA achieves good agreement with experiment. 


II. HAMILTONIAN AND METHOD OF 
CALCULATION 

The method of calculation has been described in detail 
in Ref. [50| so we give only a brief description. The 
Hamiltonian describing the NiO lattice is 

H = H 0 + Hj 

Hq = ^ ' £a diaa-diaa F H p d + H pp + Hdd 

iaa 

Hpd, = (tiajp d\ acr Pjp a + H.c.) (3) 

ict,jf3 <y 

where e.g. d\ aa creates an electron with z-spin a in the 
Ni 3d orbital a £ {xy,xz,... ,3 z 2 — r 2 } at the Ni-site 
i whereas Pjp a creates an electron in the O 2p orbital 
P £ {x,y,z} at the O-site j. The energy of the O 2p 
orbitals is the zero of energy. The terms F[ pp and F[dd 
describe hopping between two O 2p orbitals or two Ni 3d 
orbitals, respectively, and their form is self-evident. The 
parameters in this Hamiltonian have been obtained from 
a fit to an LDA band structure and are listed in Table I 
of Ref. [5(|. One noteworthy detail is that the energies 
e a have to be subject to the ‘double-counting correction’: 

—> e« — nU with U the Hubbard U - see Ref. [HJ for a 
detailed discussion. The interaction Hamiltonian FI\ has 
the form © for each Ni 3d shell, the Racah parameters 
were A = 7 eV, B = 0.13 eV C = 0.6 eV, resulting in 
the Slater integrals F° = 7.84 eV, F 2 = 10.57 eV and 
F a = 7.56 eV. An important detail is that there is a 
nonvanishing interaction only between Ni 3d orbitals in 
the same Ni ion. 

For a multiband system such as © the imaginary time 
Green’s function G(k, iw v ) and self-energy £(k, iw v ) are 
matrices of dimension 2 n or b x 2 n or b with n or b the number 
of orbitals/unit cell. In the following we will often omit 
the k, iu>v argument on these quantities for brevity. 

The starting point for the VCA is an expression for the 
Grand Canonical Potential of an interacting Fermi sys- 
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tem derived by Luttinger and Ward(5lj 


n = -i Y, eiUlu ° + (lndet(-G^ 1 )+trSG)+$[G], 

^ k,jy 

( 4 ) 


Here $[G] denotes the so-called Luttinger-Ward func¬ 
tional which was defined orginally[5l| as a sum over in¬ 
finitely many closed, connected, skeleton diagrams with 
the noninteracting Green’s function Go replaced by the 
argument of the functional, G. A nonperturbative con¬ 
struction of $[G] has been given by Potthoff }52l 153}. In 
their proof of (H]) Luttinger and Ward derived two im¬ 
portant results: first, $[G] is the generating functional 
of the self-energy 


<9$[G] 

dG a p(k, iuj v ) 


n ^/3a(k, ^^U) 


( 5 ) 


and, second, f2 is stationary under variations of S 


9So,^(k, iuj v ) 

The first of these equations can be used [47j, H3 to define 
the Legendre transform F[E] of <f>[G] via 

F[E] = $[G[E]]-±£ trSG. 

J k,j/ 


Introducing the noninteracting Green’s function Go, H 
thus can be expressed as a functional of S: 


n = -i eiUu0+ [ ln det ( _G o 1 + s ) ] + 

k,i/ 

( 7 ) 


which is known to be stationary at the exact S(w) by 
virtue of (H>|). The problem one faces in the practical ap¬ 
plication of this stationarity principle is that no explicit 
functional form of A[S] is known. 

In the framework of the VC A this problem is circum¬ 
vented as follows [4714491 ]: first, we note that <f>[G] in¬ 
volves only the interaction part Hi of the Hamiltonian 
(via the interaction lines in the skeleton diagrams) and 
the Green’s function G (via the Green’s function lines) 
- the latter, however, is the argument of the functional. 
This implies that the functional <f>[G] and its Legendre 
transform W[S] are the same for any two systems with the 
same interaction part Hi (Potthoff has derived this prop¬ 
erty without making any reference to diagrams [52t I53| ). 
In the application to NiO we accordingly consider two 
systems: System I is the original NiO lattice described 
by the Hamiltonian (O whereas System II - termed the 
reference system by Potthoff |47 h 49| - is an array of clus¬ 
ters, each of which consists of the five Ni 3d orbitals of 
one Ni ion of the original NiO lattice plus five Ligands 


or bath sites [U HH which hybridize with these. The 
single-particle Hamiltonian of such a cluster is 

Ho = ^( e d(a) dl !<T d a a + e L (a) 

a,a 

+ J2(V(a)dlJ a ^ + H.c.), ( 8 ) 

a,a 


where a £ {xy,xz.. .3 z 2 — r 2 } whereas the interaction 
part Hi for each cluster is again given by JTJ. In the Cl 
method by Fujimori and Minami the Ligand l a would be 
the linear combination of O 2p orbitals on the 6 oxygen 
ions surrounding the Ni ion under consideration which 
hybridizes with the d-orbital d a . In the case of the VGA 
the Ligands are purely mathematical objects which have 
no counterpart in the physical system and which are 
introduced solely for the purpose of constructing self¬ 
energies. Accordingly, there are no terms coupling the 
clusters centered on neighboring Ni ions in system II 
which therefore consists of completely disconnected fi¬ 
nite clusters. The crucial point is, that since the interac¬ 
tion parts of systems I and II are identical by construc¬ 
tion they have the same Luttinger-Ward functional F[S]. 
Since the individual clusters of system II are relativey 
small - they comprise 10 orbitals/spin direction - they can 
be treated by exact diagonalization and we can obtain all 
eigenstates of H — /iN within ss 20 UbT above the min¬ 
imum value. Using these the Grand Potential fl can be 
evaluated numerically (quantities with “refer to a cluster 
in the following) and the full Green’s function G(w) be 
calculated (e.g. by using the Lanczos algorithm). Next, 
G (w) can be inverted numerically for each to and the self¬ 
energy S(cc) be extracted. Thereby we have in real-space 
representation E Q = Y, a p(u) Sij where i,j are 
the indizes of the individual disconnected clusters and 
moreover Y, a p(uj) ^ 0 only if both indizes a and /? re¬ 
fer to Ni 3d orbitals. The resulting self-energy thus is 
k-independent and bears no more reference to the ficti¬ 
cious Ligands. 

Using Q and £(w) the equation © - now applied to a 
single cluster - can be reverted to obtain the numerical 
value of A[S] for the self-energy S(w). By taking the 
disgression to the reference system of clusters it is thus 
possible to generate self-energies for which the exact nu¬ 
merical value of the Luttinger-Ward functional is known. 
Next, these self-energies are used as ‘trial self-energies’ 
for the lattice i.e. we approximate 


n 


-- y 


0 iujis 0 


In det ( — G 


(-G^ + XJ) 


■ NF[±] 

( 9 ) 


where Go now is the noninteracting Green’s function of 
the physical NiO lattice and N the number of Ni-sites in 
this. 

The variation of S is performed by varying the single par¬ 
ticle parameters A, of the cluster single-particle Hamil¬ 
tonian ©, that means ed(cc), ei(a) and V(a). These 
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parameters are not determined as yet because the only 
requirement for the equality of the Luttinger-Ward func¬ 
tionals of the two systems was that the interaction parts 
H\ be identical. In this way the approximate ft © be¬ 
comes a function of the Aj, f2 = f2(Ai,... X n ) and the 
stationarity condition (jGj) is replaced by a condition on 
the A t: 



( 10 ) 


The physical interpretation would be that the VCA 
amounts to seeking the best approximation to the true 
self-energy of the NiO lattice amongst the ‘cluster rep¬ 
resentable’ ones. Since its invention by Potthoff the 
VCA has been applied to study the Hubbard model in 
various dimensions ImMchI models for 3d transiton metal 
compounds and interacting Bosons [70lf71j. 

For the present application to NiO and in the paramag¬ 
netic case cubic symmetry reduces the number of param¬ 
eters A,; to be varied to only six: for each a the Hamil¬ 
tonian © contains 3 parameters and there is one such 
set for the e g orbitals and one for the t 2 g orbitals. The 
equation s yste m (fTU|) is solved by the Newton method, 
see Refs. I5C |, [68j for details. The proposal of Balzer 
and PotthofT[5j| to use rotated and rescaled coordiante 
axis for the calculation of the derivatives of U with re¬ 
spect to the A i turned out to be of crucial importance for 
successful Newton iterations. 

The paramagnetic phase of NiO was studied in some de¬ 
tail in the preceeding paper Ref. [50]. There are only a 
few differences as compared to this study: first, a reduced 
value of the Racah-parameter A = 7.0 eV (A = 8.25 eV 
in Ref. 0) and consequently a readjustment of the Ni 
3d-orbital energy e<z = — 52 eV (whereas ed = — 62 eV in 
Ref. [50j ) to account for the different double-counting 
correction. Moreover, the Ni 3d-to-0 2p hopping param¬ 
eter (pda) was increased by 10% to —1.4178 eV. 

Due to improved computer power it was now moreover 
possible to optimize all relevant single particle parame¬ 
ters of the octahedral cluster. Thereby it turned out that 
V(t 2 g) = 0 is a stationary point irrespective of the val¬ 
ues of the other parameters. The energy of the t 2 g -Ake 
Ligand, eL(t- 2 g) then is irrelvant so that only the four 
parameters ed(e g ), ei(e g ), V(e g ) and Cd(t 2 g ) remain to 
be solved for. The resulting paramagnetic solution, how¬ 
ever, is stationary with respect to all six possible param¬ 
eters. Some results for the paramagnetic phase will be 
presented later in comparison to the antiferromagnetic 
one. 


III. MAGNETIC SUSCEPTIBILITIES AND 
ANTIFERROMAGNETIC TRANSITION 


We discuss the staggered and uniform magnetic suscep¬ 
tibility. Within the VCA the Grand Canonical Potential 
may be thought of as being expressed as a function of a 


number of parameters 

0 = 0(Ci,... Cm, Ai,..., A n ), (11) 

were the Q are the parameters of the physical lattice 
system - such as the physical hoppig integrals and orbital 
energies or certain external fields - and the A i are the 
single-electron parameters of the reference system which 
parameterize the self-energy. We assume that amongst 
the Q there is also a uniform or staggered magnetic field h 
along the 2 -direction. This implies that the values of all 
single-particle parameters of the reference system must 
be taken as spin-dependent: 

A^tr — Aq_|_ -(- sign(a)Ai } _, 

which results in a spin-dependent self-energy, E-j-(w) ^ 
£_j_(w). For a staggered field we switch to the antiferro¬ 
magnetic unit cell and assume that the A^._ have opposite 
sign at the two Ni-ions in this cell. This means that the 
self-energy for an f-electron is £-|-(w) at the first Ni-ion 
and Ej.(w) at the second Ni-ion in the antiferromagnetic 
cell and vice versa for a ^.-electron (k-sums now have to 
be performed over the antiferromagnetic Brillouin zone). 
For a uniform field we retain the orginal unit cell and use 
the spin dependent self-energy at the single Ni-ion in this 
cell. 

We assume that we have found a stationary point A* for 
h = 0 i.e.: 


on 

dX i 



( 12 ) 


for all i and denote the Grand Potential for this solution 
by D 0 . Since this is the paramagnetic stationary point 
all spin-odd parameters A^- are zero. Upon applying a 
small finite h in the lattice system , f1 therefore can be 
expanded as 


fi — Do + — y ] Xi A, 


i i3 A? ~b 'y ] A i Bi h + — C 


(13) 


where the shifts Xi = Xi — A* and A, B and C are 2 nd 
derivatives of U at the point Aj = A*, h = 0. There are no 
terms linear in the A i because of (11211 and there is no term 
linear in h because U must be an even function of h. The 
2 nd derivatives can be evaluated numerically whereby the 
fact that the Xi are parameters of the reference system 
whereas the staggered field h is one of the Q in (fTTh causes 
no problem. Moreover, all derivatives are to be evaluated 
in the paramagnetic phase, so no calculation in a finite 
field is necessary. Demanding stationarity we obtain for 
the shifts Xp. 

OQ 

Tpr = Ai,jXj + B t h = 0, (14) 

O Aj 

L n 
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FIG. 1: Staggered (top) and uniform (bottom) susceptibility 
of NiO. The symbols give the calculated values, the lines are 
the expressions m and inl¬ 


and reinserting into (1131) we obtain fl as a function of 

h \ n 

n(h) = n 0 - x 

X = E - c - (is) 


f 2 must be invariant under a simultaneous sign change of 
h and all spin-odd parameters A i- so that A^+ = 0 for 
all i and the sums over i and j in m extend only over 
the spin-odd parameters. 

Figure [l] shows the staggered and uniform suscepti- 
bilites obtained in this way as a function of temperature. 
Thereby the term 


was added to the lattice Hamiltonian, so that the physical 
susceptibilities are given by 


2 d 2 n 

XQ ~ “Ms 

The staggered susceptibility can be fitted accurately by 


Xs(T) 


Ml C s 
T — T n 


(16) 


with Tjy = 481.03 K and C s = 24805 K ■ eV 1 whereas 
the uniform susceptibility can be fitted by 


Xo(T) 


m|Cq 
T + Qgw 


(17) 


with Sew = 491.62 K and C 0 = 30923 K • eH" 1 . The 
divergence of Xs(T) at T )v is due to an eigenvalue of the 
Hessian A crossing zero at Tn- 

Interestingly the Curie-Weiss temperature Qcw is some¬ 
what higher than the Neel temperature T'n- This can be 
understood as a consequence of the type-II antiferromag¬ 
netic structure and (weak) antiferromagnetic exchange 
between Ni-ions connected by a 90°-degree Ni-O-Ni bond 
(an example would be two Ni-ions at distance (a, a, 0)). 
For such a pair of Ni ions there is a competition between 
the direct antiferromagnetic exchange (mediated by the 
direct Ni-Ni hopping as described by Slater-Koster pa¬ 
rameters such as ( dda )) and the ferromagnetic exchange 
due to Hund’s rule coupling on oxygen. Let us assume 
that the net exchange constant between such a pair of 
Ni ions is antiferromagnetic (this is certainly true for the 
present calculation which does not include Hund’s rule 
exchange on oxygen). Then, any given Ni ion has 12 
neighbors of that type and in the type-II antiferromag¬ 
netic structure the ordered moment of one half of these 
neighbors is parallel to the moment of the ion at the 
center whereas it is antiparallel for the other half. The 
exchange fields due to these 12 neighbors therefore cancel 
and the Neel temperature is determined solely by the an¬ 
tiferromagnetic superexchange with the 6 neighbors con¬ 
nected by 180°-degree Ni-O-Ni bonds. On the other hand 
the parameter Qcw is a measure as to how strongly the 
antiferromagnetic exchange between spins opposes a uni¬ 
form ferromagnetic polarization. For the case of a uni¬ 
form ferromagnetic polarisation, however, the exchange 
fields from all 12 (a, a, 0)-like neighbors are parallel and 
therefore do contribute to Qcw- The discrepancy be¬ 
tween Tn and O obtained in the VCA calculation thus 
is to be expected. The experimental value Q exp = 2000 
K was given in Ref. [72] but being almost four times the 
experimental Neel temperature this appears somewhat 
high. 

For a Heisenberg antiferromagnet the constant C would 
be given by 


H m = -h E e lQ Ri (n it - n i4 .) 

i 


c = 


S(S + 1) 
3 ks 


CgSeff ) 2 
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Using the values from the fits in Figure |T] with S = 1 
we obtain the reasonable values S e ff = 1.000 for Xo and 
S e ff = 0.895 from \s- The smaller value for S e ff in the 
staggered case likely is due to the fact that in the com¬ 
putation of xo the magnetic field is applied to all orbitals 
in the unit cell whereas it acts only on the Ni 3d orbitals 
in the case of XV 

Lastly we point out an interesting feature of these re¬ 
sults: the matrix A of second derivatives of Q with re¬ 
spect to the spin-odd parameters Aq_ in (fl5l) obviously 
is the same in the case of staggered and uniform suscep¬ 
tibility. Only the quantities B and C are different. Still, 
the resulting susceptibilities have a completely different 
but physically reasonable temperature dependence. 


IV. ANTIFERROMAGNETIC PHASE 

To obtain the results presented so far only the para¬ 
magnetic solution was needed. If we want to discuss the 
antiferromagnetic phase itself we need to find station¬ 
ary points with spin dependent parameters A i and this 
doubling of the number of parameters complicates the 
numerical problem of finding the stationary point. We 
recall that in the parmagnetic case a total of 4 parame¬ 
ters were varied: ed(e g ), eL(e g ), V(e g ) and ed{t 2 g ) (more¬ 
over, V(t 2 g ) = 0 always was a stationary point and with 
that value the last parameter e.L(f 2g ) is irrelevant). Dou¬ 
bling all of the nonvanishing parameters would result in 
a total of 8 parameters. Using the Newton method this 
is still numerically manageable but it turned out that 
a more severe problem appears. While it might seem 
that the more parameters \ one is varying the better 
an approximation for the self-energy results, calculations 
showed that the opposite is true. Introducing too many 
symmetry-breaking parameters Aleads to unphysical 
solutions - one example is discussed in detail in Appendix 
A. In fact it turned out that retaining more than 2 - out 
of the 4 possible - Aleads to unphysical solutions. Ac¬ 
cordingly, in the following we present solutions obtained 
with 6 parameters A i. In choosing the Ato be kept 
we heuristically use the staggered susceptibility Xs as a 
guidance. Namely we can restrict the set of the Aq_ in 
the expression (1151) for \s to only 2 and examine which 
combination still gives a x s which is closest to the one 
obtained with the full set of 4 Aq_. It turned out that 
retaining only the spin-odd part of the e g -like hopping in¬ 
tegral V_(e g ) and the f 2g -like d-level energy e< 2 ,-(i 2g ) the 
staggered susceptibility Xs - and in particular the Neel 
temperature - remain practically unchanged. This ap¬ 
pears plausible because the physical mechanism that sta¬ 
bilizes antiferromagnetism in NiO is the enhanced hop¬ 
ping for the e g electrons of one spin direction along the 
180°-degree Ni-O-Ni bonds connecting sites on different 
sublattices. A spin dependent d-level-to-Ligand hopping 
then clearly is the best way to simulate this effect in 
a cluster with only a single Ni-ion. Since we have set 
U(f 2g ) = 0 in the paramagnetic case, £d(t 2g ) moreover 



Temperature [K] 

FIG. 2: Grand Potential per Ni ion for the various solutions 
described in the text as a function of temperature. 

Do = —356.042 eV is an arbitrary reference energy. 


is the only remaining parameter pertinent to the t 2g or¬ 
bitals. In the following, the resulting solution will be 
referred to as AF-I. 

In addition there is a second type of antiferromagnetic 
solution where V+(e g ) = U_(e g ) so that the hopping for 
one spin direction of e g electrons is exactly zero. In this 
case the only remaining spin-odd parameters to be var¬ 
ied are £d,-{t 2 g ) and ed,-(e g ) (since one spin direction 
of the e g electrons has zero hopping, the spin splitting 
of the e g -like ligand, ez /j _(e g ) is irrelevant and can be 
set to zero). It should be stressed that in this case H 
is stationary also with respect to V_(e g ), and the values 
V+(f 2g ) = V_(f 2g ) = 0 remain stationary as in the para¬ 
magnetic case. Despite the fact that only 6 parameters 
are actually varied, the corresponding solution therefore 
is stationary with respect to all 12 possible parameters 
of the cluster. This solution will be referred to as AF-II. 

Figure [2] shows H as a function of temperature for the 
paramagnetic, AF-I and AF-II solutions. For the para¬ 
magnetic phase to good approximation U(T) = Do — 
ksTlog{ 3) where kslog{ 3) is the entropy due to the 3-fold 
degenerate 3 A 2g ground state of a single Ni 2+ ion with 
configuration t 2g e g . At Tjv the solution AF-I branches off 
as would be expected for a 2 nd order phase transition. At 
237.5 K there is a crossing with finite difference of slopes 
between the AF-I and the AF-II solutions. This would 
imply a 1 st order phase transition which most probably 
is unphysical. Rather, this may be the way in which the 
VC A approximates a continuous but rapid change of the 
electronic state. 

We discuss the phase transition at Tm- As usual we in¬ 
troduce the staggered magnetization m s 


m s = — 


dfl 

dh s 


(18) 
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FIG. 3: Staggered field versus staggered magnetization below 
and above Tn- The lines are the cubic fits discussed in the 
text. 


and switch to the Legendre-transformed potential 
f l'(m s ) = Cl(h s ) + m s ■ h s which - using (1151) - is 

fi 'K) = \A Xj 1 rn 2 s + 0(m 4 ). 

In the present situation where the chemical potential is 
within the sizeable insulating gap this equals the Gibbs 
free energy up to an additive constant. Comparison with 
(1161) shows that the lowest order term in the expansion 
of f Y(m s ) has the form expected from Landau theory: 

G'(m s ) = a(T - T N )m 2 + ^ m 4 
with a = 1/(2 C s ). To carry this further we use 

h s = —— : 2 a(T — T N )m s + 2 bm 3 . 
am s 

The dependence of h s on m s can be obtained within the 
VCA by increasing the staggered field h s in small steps 
starting from h s = 0, thereby always using the converged 
AF-I solution for the preceding step as starting point for 
the Newton algorithm for the next h s . After convergence, 
m s is obtained from m- The resulting curves are shown 
in Figure [3] At 470 K the magnetization is opposite to 
the field because below T)v the paramagnetic state and 
hence also the states ‘close’ to it are unstable and in fact 
this behaviour is precisely what is expected from Landau 
theory. Namely h s (m s ) can be fitted by a third order 
polynomial h s = c\m s + C2TO^ with ci = —0.484 meV, 
C2 = 2.339 rneV at 470 K and c\ = 0.322 meV, C2 = 
2.458 rneV at 490 K. From the rather similar values of 
C2 we can conclude that b ss 1.20 meV. Evaluating c\ 
from ci = 2a(T — Tn) gives ci = —0.447 meV at 470 
K and ci = 0.362 meV at 490 K, reasonably consistent 
with the values extracted from m s (h s ). Figure U shows 




FIG. 4: Top: Difference of the Grand Potential in the param¬ 
agnetic phase and in the antiferromagnetic phase versus tem¬ 
perature (dots) and fit (line). 

Bottom: Ordered moment versus temperature (dots) and fit 
(line). 

the difference G paT , Q — Haf ■ Since both are calculated in 
zero external field we have O' = fl and we expect 

n para (T)-n A F(.T) = A(T — T n ) 2 - 

The fit gives A = 1.746 10~ 4 meV K~ 2 whereas us¬ 
ing a 2 /{2b) = 1.694 10~ 4 meV K~ 2 . Figure []] also 
shows the ordered moment m s = (rid, f) — ( n d,i) ver¬ 
sus temperature. Close to Tn this can be fitted by 
m s (T) = B^JTn - T with B = 0.131 K~ 4 ! 2 (for com¬ 
parison: ( a/b) = 0.130 I\ 4 / 2 ). The symmetry¬ 

breaking parameters V-(e g ) and ed,-(e g ) have a similar 
T-dependence oc \/Tn — T close to Tn- All in all the 
phase transition is described well by Landau theory and 
the VCA allows to extract the parameters of the theory 
from the original Hamiltonian. 

Next we consider the specific heat C(T) which is shown 
in Figure [S] The top panel shows experimental data for 













FIG. 5: Top: Experimental spefific heat of NiO as obtained 
by different authors [734751 and phonon contributionf76|. 
Bottom: Experimental specific heat after subtraction of the 
phonon contribution compared to the VCA result. The ver¬ 
tical line gives the temperature where the solutions AF-I and 
AF-II cross, see Figure [2] 


NiO taken from Refs. j73l475| . Also shown is the phonon 
contribution which was calculated from the phonon spec¬ 
trum measured by inelastic neutron scattering at room 
temperature [l^]. The electronic heat capacity C e i is com¬ 
pared to the VCA result in the lower part of the figure. 
The VCA of course does not reproduce the divergence 
of C e i at the ordering transition which follows the crit¬ 
ical exponents for a 3D Heisenberg antiferromagnetj75j. 
Apart from that the result from the VCA is roughly con¬ 
sistent with the measured values, especially the solution 
AF-II agrees quite well with the data at low temperature. 
In fact, this solution appears to match the experimental 
data considerably better up to ~ 0.7 T'n- This may be 
an indication that in NiO this solution is realized even at 
temperatures well above the transition between AF-I and 
AF-II at 237.5 K (which is indicated by the vertical line 
in Figure [5j). In fact, as will be seen below, the single 
particle spectral function for this solution matches the 
experimental photoelectron spectra - which are usually 
taken at room temperature - quite well. 


The apparently large value of C e i above Tjv probably is 
an artefact: the phonon spectrum was measured at room 
temperature that means where the lattice was deformed 
by magnetostriction. Above T/v this deformation is ab¬ 
sent and the phonon spectrum may change so that using 
the low temperature phonon spectrum gives an incor¬ 
rect estimate for the phonon contribution. The electronic 
heat capacity per mole obeys the sum rule 

J q h — 7 ^—^ dT = i?log(3) 

where Th is well above the Neel temperature. It turns 
out that the experimental data exhaust this sum rule at 
T h « 660 K. 

Finally, Table U compares the parameter values for the 



AF 

Para A 

e d(£g) t 

-51.4113 

-51.5278 

C-di&g') 4 r 

-50.9896 

-51.5278 

V(e g ) t 

-1.4880 

-1.7260 

V(e g ) ( 

0.0000 

-1.7260 

ez,(e 9 ) t 

4.3633 

3.4089 

ez,(e g ) 4- 

4.3633 

3.4089 

ed(t2g) t 

-51.1777 

-51.2256 

ed(t2g) ( 

-51.0794 

-51.2256 

12 + fiN 

-244.0363 

-244.0235 -0.0128 

( H ) 

-244.0344 

-244.0046 -0.0299 

S/k B 

0.1090 

1.0977 -0.9888 

(Ho) 

-457.9924 

-457.3822 -0.6102 

(Hi) 

213.9579 

213.3776 0.5803 

(H P d) 

-3.5283 

-3.4380 -0.0903 

(Hpp) 

-0.1915 

-0.1757 -0.0158 

(H dd ) 

-0.0030 

-0.0034 0.0004 

(«e 9 ,t> 

0.2247 

1.0945 -0.8698 

( n e g ,l) 

1.9738 

1.0945 0.8793 

( n d) 

8.1985 

8.1890 0.0095 

(nt 2g ,t) 

3.0000 

3.0000 0.0000 

(nt 2a ,i) 

3.0000 

3.0000 0.0000 


TABLE I: Top: Parameters of the AF-II (AF) and param¬ 
agnetic (Para) solution at 200 Kelvin. A=AF-Para. Center 
part: Various contributions to 12 (energies in eV per Ni) Bot¬ 
tom part: Occupation numbers of the Ni 3d shell. 


paramagnetic and AF-II solutions at 200 Kelvin, as well 
as various observables. The expectation value of any 
term H part in the Hamiltonian can be calculated by re¬ 
placing H part —I (Hpart and using 


(Hp art ) = I? 


C =1 


The numerical calculation of the derivative thereby is 
simplified considerably by taking into account that due 
to the stationarity condition for the At their variation 
with ( can be neglected [56j. The antiferromagnetic unit 
cell and k-mesh were used also for the paramagnetic cal¬ 
culation to avoid artefacts. 

A somewhat surpising feature is that the quite different 














9 


cluster parameters for the two different solutions give 
only slightly different f2. As expected, f 1 is lower for 
the antiferromagnetic phase due to the lower energy ( H). 
This is partly compensated by the almost complete loss 
of entropy in the AF phase but at the low temperature 
considered this does not result in a higher fl. As already 
mntioned the entropy in the paramagnetic phase is close 
to S/ks = log( 3) = 1.0986 as expected for a system of 
localized S = 1 spins. While this may seem trivial it 
should be noted that the spin degeneracy can only be 
reproduced if a spin rotation invariant Hamiltonian is 
used. Discarding the off-diagonal matrix elements of the 
Coulomb interaction m breaks the spin-rotation symme¬ 
try so that the entropy of the paramagnetic phase cannot 
be obtained correctly. 

The considerably lower value of (Ho) in the AF phase 
comes about because the d-occupancy increases slightly, 
by 0.0095. With the orbital energy = —52 eV this low¬ 
ers (Ho) by —0.494 eV. This is more than compensated, 
however, by the increase of the Coulomb energy (Hi) in 
the AF phase by 0.5803 eV. Eventually, the energy in 
the antiferromagnetic phase becomes lower due to the 
dp-hybridization which is lowered by —0.0903 eV in the 
AF phase. Consistent with the theory of superexchange 
the driving force behind the antiferromagnetic ordering 
is not the lowering of the Coulomb energy but a gain of 
kinetic energy. Interestingly there is also a significant - 
on the scale of the change of ( H ) - gain in the direct O 
2p-0 2p hopping energy ( H pp ) in the AF phase. This 
is due to the increased charge transfer to Ni 3d which 
reduces the Filing of the O 2p orbitals and thus allows 
for enhanced O 2p-0 2p hopping. This contribution is 
missed in models for superexchange which consider only 
a single Ni-O-Ni bond. 

To summarize the results of the two preceding sections: 
the description of the magnetic properties and phase 
transition of NiO as obtained by the VCA is very simi¬ 
lar to what would be obtained from a simple mean-field 
treatment of a localized spin system. It should be noted, 
however, that the Hamiltonian does not contain any ex¬ 
change terms nor is there any molecular field in the phys¬ 
ical system. Rather, the Hamiltonian is the one for the 
full NiO lattice ([3]) , the self-energy is calculated in a clus¬ 
ter containing a single Ni-ion and the coupling between 
Ni-ions is solely due to the lattice kinetic energy. Still, 
the VCA even captures subtle details such as the pres¬ 
ence of different exchange channels as manifested by the 
different values of Tjv and Qcw- It should also be kept in 
mind that all parameters in the original Hamiltonian are 
several orders of magnitude larger than the differences 
in energy in Table Q] but still the calculated Neel tem¬ 
perature is quite close to the experimental value. The 
VCA thus appears successfid in correctly extracting the 
low energy scales relevant for ordering phenomena and 
thermodynamics from the high energy scales of Hubbard 
U, charge transfer energy, and hopping parameters. 



FIG. 6: fc-integrated spectral function A(lo) for the antifer¬ 
romagnetic solution AF-II at 200 Kelvin compared to exper¬ 
imental spectra. 


V. PHOTOELECTRON SPECTRA IN THE 
ANTIFERROMAGNETIC PHASE 

We proceed to a discussion of ‘high energy physics’ 
and consider the single particle spectral function. Figure 
[6] shows the fc-integrated spectral density 

A(w) = G ata (k,u]+ irj), (19) 

hoc 

where the sum over a runs over either the Ni 3d or the O 
2p orbitals. This was calculated for the antiferromagnetic 
solution AF-II at 200 Kelvin. Experiments are usually 
done at room temperature but the specific heat data in 
Figure [5] suggest that the AF-II solution may be relevant 
also at higher temperature. The Figure also shows dif¬ 
ferent experimental angle-integrated spectra: First, hard 
x-ray photoelectron spectro scop y (HAXPES) with a pho¬ 
ton energy of hv = 6500 eV[73] - at this energy the pho¬ 
toionization cross section for Ni 3d is approximately 10 
times larger than that for O 2p so that predominantly the 
Ni 3d density of states is observed. Second, an x-ray pho¬ 
toemission (XPS) spectrum taken with a photon energy 
of hv = 67 eV[ll]. This is close to the Ni 3p—> 3d absorp¬ 
tion threshold so that the satellite region around —10 eV 
is resonantly enhanced[79l|. And, third, an x-ray emis¬ 
sion (XES) spectrum which shows predominantly the O 
2 p density of states (80], 

Several peaks in the theoretical spectra can be identified 
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in the various experimental spectra: these are the peaks 
A and B which have Ni 3d character and thus appear 
in HAXPES and XPS (although at hv = 67eV the peak 
A is anti-resonantly suppressed). The peak A also has 
some oxygen admixture so that together with peak C it 
can also be seen in the XES spectrum. Peak C also cor¬ 
responds to a weak feature observed in HAXPES and in 
the XPS spectrum. Peak D can be seen both in HAX¬ 
PES and XPS and the tail of the XES spectrum towards 
negative energy also shows an indication of the additional 
shoulder E which corresponds to a similar feature in the 
theoretical O 2p spectrum. Finally the rather broad fea¬ 
ture F can be seen very well in the XPS spectrum. By 
and large there is good agreement between calculated and 
measured spectra. It has to be kept in mind, however, 
that as far as the Ni 3d density is concerned a similar 
degree of agreement has been obtained earlier by Fuji¬ 
mori and Mina,mi [ 20 I ] and van Elp et aZ.[26j by the clus¬ 
ter method. The discussion so far shows mainly that as 
far as angle-integrated spectra are concerned the VGA 
‘inherits’ the accuracy of the cluster method. 

We therefore turn to the quantity which allows for the 
most detailed comparison to experiment, namely the k- 
resolved spectral density 

A(k,w) = — 'y^G a>a (k,u + irj), (20) 

7r z ' 

a 

where the sum over a now runs over both, the Ni 3d 
and the O 2p orbitals. The dispersion of peaks in A(k, ui) 
can be compared to the band structure as measured in 
angle resolved photoemission spectroscopy (ARPES). To 
date there are two ARPES studies of NiO, one by Shen et 
al. [11] and the other by Kuhlenbeck et al. |82j. Shen et al. 
give three different sets of data points: the bands from T 
to X = (221,0,0) (i.e. the (1, 0,0)-direction) measured in 
normal and off-normal emission and the bands from T to 
Xi = (qf,qf,0) (i.e. the (1,1,0)-direction)) measured 
in off-normal emission. The two data sets along P — X 
agree for some bands but differ for others due to ma¬ 
trix element effects. If a given band is observed in any 
one experimental geometry it obviously does exist and 
if it is not observed in another geometry this can only 
be be a matrix-element effect. The true band structure 
along T — X thus should comprise at least the superpo¬ 
sition of the two sets of bands for normal and off-normal 
emission. Figure [7] compares A(k,u) and the respective 
experimental band dispersions. For both directions the 
top of the band structure is formed by a complex of sev¬ 
eral closely spaced bands with high spectral weight in 
the range —0.5 —> —1.5 eV, labeled a in the Figure. In 
the angle-integrated spectrum in Figure [G] these bands 
produce the intense peak A. The high spectral weight of 
these bands can also be seen in the experimental spectra 
in Figures 7 and 8 of Shen et al.. Below this group of 
bands there is a gap of approximately 1 eV. In the range 
—2.5 —> —3.5 eV there are several essentially dispersion¬ 
less bands with weak intensity, labeled b. In the angle 
integrated spectrum in Figure |B] these bands produce the 



FIG. 7: Single particle spectral density A(k, w) for the antifer- 
romagneric solution AF-II at 200 Kelvin compared to ARPES 
data by Shen et al.[8l|. The three panels show the data for 
off-normal emission along F — X (top), normal emission along 
r — X (middle) and off-normal emission along F — A'i (bot¬ 
tom). 


weak feature B. Shen et al. resolved two such bands 
along (1,0, 0) but only one along (1,1,0) - since the dis¬ 
persions must match at T there probably are more than 
one of these dispersionless bands also along (1,1,0). 
Such (nearly) dispersionless bands can be seen also at 
even more negative energies, but there they are super¬ 
imposed over and mix with the strongly dispersive O 
2p derived bands which results in more fc-dependence. 
First, there is the dispersive band c along (1,0,0) and e 
along (1,1,0). In experiment band c shows a relatively 
strong upward bend near X - this may indicate that there 
rather a part of the dispersionless band which starts at 
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FIG. 8: Single particle spectral density A(k, u) for the anti- 
ferromagneric solution AF-II at 200 Kelvin compared to the 
ARPES data by Kuhlenbeck et al. f&|. 


r at « —4.6 eV has been observed. This dispersionless 
band and its ‘avoided crossing’ with a strongly disper¬ 
sive O 2p derived band c as predicted by the VGA may 
also have been observed near T, see the region labeled / 
in the top panel. Along (1,1,0) this dispersionless band 
can be followed over the full k-range, see the band la¬ 
beled d in the bottom panel. In normal emission (middle 
panel) it moreover becomes apparent that the strongly 
downward dispersing O 2p band indeed splits into two 
bands with opposite curvature - see the region labeled g 

- which would be similar to the VCA bands. Lastly, at 
« —6.6 eV another dispersionless band - labeled g and h 
in the middle and bottom panel - is observed which would 
correspond to the nearly dispersionless band which starts 
out from r at slightly below —6 eV (and which gives rise 
to the peak D in the angle-integrated spectrum in Figure 

Some of the above interpretations are corroborated in 
Figure |5] which shows a comparison to the band structure 
deduced by Kuhlenbeck et al. along T — There 

the dispersionless bands a, b and c obviously correspond 
to the bands with the same labels in the data by Shen et 
al. - and the corresponding bands predicted by the VCA. 
Particularly interesting is the band portion labeled d in 
Figure[S]which also shows a peculiar downward curvature 
and corresponds exactly to the part labeled g in Figure [7] 

- which in turn had some counterpart in the VCA bands 
(plotting the two experimental band structures on top 
of each other shows the exact correspondence of these 
two bands). Finally, the dispersionless band portion e is 
precisely the continuation towards T of the dispersionless 
band labeled h in Figure Q]- which also has its counter¬ 
part in the VCA band strcuture. 

Generally speaking for all bands which should be easy 
to observe because they either have a high intensity or 
are relatively isolated from other bands there is an es¬ 
sentially one-to-one correspondence between VCA and 
experiment. The VCA predicts a multitude of disper¬ 
sive low-intensity bands below ss —4 eV and only a few 



FIG. 9: k-integrated spectral density for the antiferromag¬ 
netic AF-II solution at 200 K and the paramagnetic solution 
at 520 K. Also shown are experimental spectra by Tjernberg 
et al. measured below and above the experimental Neel tem¬ 
perature of 523 K. 


of these seem to have been observed. Combining the 
three experimental spectra along T — A' indicates, how¬ 
ever, that the experimental band structure in this en¬ 
ergy range does not consist just of the three strongly dis¬ 
persive O 2p-derived bands obtained by band structure 
calculations [2| which are also predicted by many ‘cor¬ 
related’ calculations as well [Til fl^ |. Rather, additional 
bands, both dispersive and non-dispersive, appear to be 
observed. In the next section it will be shown that the 
dispersionless bands are in fact the very fingerprints of 
the atomic multiplets in the ARPES spectra. 

To conclude this section we discuss the temperature de¬ 
pendence of the spectra. The bottom part of Figure [9] 
compares the angle integrated Ni 3d-like spectral func¬ 
tion for the AF-II solution at 200 K and for the para¬ 
magnetic solution at 520 K (the paramagnetic solution 
has practically no temperature depedence). While the 
spectra have very similar overall shape there are small 
differences. The gap between the large peaks A and D 
is filled with weight and the relatively well-defined peak 
B more or less disappears in the paramagnetic phase. A 
rather strong redistribution of weight occurs in the satel¬ 
lite region were weight disappears between —11 eV and 
—8 eV and a new strong peak grows at « —8 eV. 
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Figure [9] also shows the experimental spectrum taken by 
Tjernberg et aZ. (T^] at 615 K - which is well above Tjv - 
and a modified version of the spectrum at 300 K. More 
precisely, the 300 K spectrum was broadened by convo¬ 
lution with a Gaussian to simulate the enhanced thermal 
broadening and fitted to the high-temperature spectrum 
whereby both spectra were normalized to unity [781]. Due 
to the relatively low photon energy hv = 65 eV the ex¬ 
perimental spectrum also contains a considerable amount 
of O 2 p weight, which gives rise to the intense peak C 
(compare Figure [ 6 ]). Accordingly Figure [HI also shows the 
sum of Ni 3d-like and O 2p-like spectral functions. In the 
experimental spectrum the peak A looses weight in the 
paramagnetic phase whereas the opposite is predicted by 
the calculation. In experiment the spectral weight in the 
energy range between the two large peaks A and C in¬ 
creases in the paramagnetic phase and a similar change 
occurs in the theoretical spectra were the relatively well- 
defined gap between the peaks A and C is partly filled 
in the paramagnetic spectrum, although the effect seems 
less pronounced in experiment.In experiment the peak C 
looses a small amount of spectral weight and is shifted to 
slightly less negative energy in the paramagnetic phase. 
A similar tendency can be seen in the theoretical spec¬ 
tra but considerably exaggerated. Lastly, in experiment 
the spectral weight increases slightly at various positions 
in the satellite region below —8 eV but no decrease is 
observed anywhere. In contrast to this in the theoreti¬ 
cal spectra there is a drastic change in the satellite region 
where a considerable amount of weight disappears around 
— 10 eV and a new strong peak appears at approximately 
—8 eV. Summarizing, the VCA is only partly success¬ 
ful in predicting the changes of the photoemission spec¬ 
trum across the Neel temperature. It has to be kept in 
mind, however, that the photon energy of 65 eV used 
in the experiment is close to the 3 p —» 3d absorption 
threshold so that the satellite (peak A) are resonantly 
enhanced (antiresonantly suppressed). In fact, the in¬ 
tensities of the various peaks are quite different from 
the HAXPES spectrum in Figure [ 6 ] Accordingly, ad¬ 
ditional effects may come into play which determine the 
intensity of these features and this might be one expla¬ 
nation why discrepancies with theory occur precisely for 
peak A and the satellite. In any way some of the ob¬ 
served changes with temperature - or actually: between 
antiferromagnetic and paramagnetic phase - appear to 
be reproduced qualitatively by the VCA. Lastly, Figure 
[9] also shows a somewhat surprising difference between 
the single particle spectra in the antiferromagnetic and 
paramagnetic phase: namely the insulating gap in the 
paramagnetic phase is larger than in the antiferromag¬ 
netic phase. More precisely, the peak-to-peak distances 
are 4.65 eV and 4.05 eV so that the gap increases by 
« 10 %. So far the temperature dependence of the insu¬ 
lating gap in NiO has not been studied experimentally. 
As will be discussed in the next section, however, there is 
a clear physical reason reason for this discrepancy namely 
the fact that the mechanism which opens the insulating 



to [arb. u.] 

FIG. 10: Graphical solution of equation (1221) . The self-energy 
has a single pole at £ = 1 . 

gap in the two phases is quite different. 

VI. DISCUSSION OF THE SELF ENERGY 

We discuss some of the results presented in the pro¬ 
ceeding section from the ‘self-energy perspective’. Lut- 
tinger has shown[83j] that for a single band system the 
self-energy has a spectral respresentation of the form 

E(u;) = ?7 + ^ ( 21 ) 

i u & 

where 77 , eq > 0 and Q are real. In the following we as¬ 
sume these parameters to be k- independent. The equa¬ 
tion for the poles of the Green’s function reads 

uJ k - e k = V —(22) 

l Wfc - C i 

where for brevity of notation we have replaced efc +77 — 
/i —> efc. Let us first assume that we have only a single 
pole with a large weight, Figure flOl shows the resulting 
£(w) for real ui. Also shown is the noninteracting density 
of states for the band e*,, the two straight lines correspond 
to oj — e_ and u> — e + where e_ and e + are the bottom 
and top of the noninteracting band ek- The intersections 
of these lines with E(w) give the solutions of the equa¬ 
tion (|22D and there is one solution for any k in between 
these. A single isolated pole of E(w) thus splits the non¬ 
interacting band into the two Hubbard bands and opens 
a gap in the spectral function. Such an isolated pole 
with a residuum oc N° in the self-energy obviously is the 
very essence of a Mott insulator - this can also be seen in 
the self-energy of the 2-dimensional Hubbard model [S4j. 
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FIG. 11: Top: e g -like self-energy E(cu) and k-integrated d-like 
spectral function A(uj) for the paramagnetic solution at 520 
K. Bottom: same for the AF-II solution at 200 K. The two 
different self-energies for the antiferromagnetic case refer to 
the two spin-directions. The scale for the self-energy is the 
same for both panels. 

Figure fill shows the e ff -like self-energy for the paramag¬ 
netic solution at 520 K and for the antiferromagnetic AF- 
II solution at 200 K as well as the respective k-integrated 
d-like spectral function A(w). For the paramagnetic so¬ 
lution there is indeed an isolated intense peak of the self¬ 
energy - labelled G in the Figure - within the insulating 
gap. In contrast no such ‘gap-opening-peak’ exists in the 
self-energy for the antiferromagnetic solution. There the 
mechanism which opens the gap is the different values of 
the additive constants 77 -j. and ij 4 . in the spin-dependent 
self-energy, which has the effect of an oscillating potential 

which opens a gap in the same way as in spin-density- 
wave mean-field theory. For the AF-II solution at 200 
K ?7t( e s) = 59.89 eV 77 i(e g ) = 50.69 eV (whereby the 
large average (rj^ + rji )/2 ss 55.3 eV cancels the double¬ 
counting correction to the Ni 3d-level energy). Since the 
insulating gaps in the paramagnetic and antiferromag- 




Energy [eV] 



Energy [eV] 


FIG. 12: Temperature variation of Vsdw = (??t — wi)/2 and 
the e g -like self-energy for the two spin directions in the AF-I 
solution. 


netic phase are created by different mechanisms it is not 
too surprising that they have different values (see Fig¬ 
ure [9]). For completeness we mention that the t 2 g -like 
self-energy has no such ‘gap-opening-peak’ in either the 
paramagnetic or insulating phase. This is special for NiO 
because in the 3 A. 2 9 ground state of d 8 in cubic symme¬ 
try, which is tt g e % the f 2g orbitals are completely filled 
and thus comprise a ‘band-insulating’ subsystem. Very 
probably this is also the reason why V{t- 2 g ) = 0 is a sta¬ 
tionary point. 

The AF-I solution which branches off the paramagnetic 
solution at Tjv and crosses with the AF-II solution at 
237.5 K (see Figure [2jl interpolates between these two 
types of insulating gap: Figure [T5] shows £(w) for this 
solution for different temperatures. With decreasing tern- 
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to [arb. u.] 

FIG. 13: Top: Spectral density of the ‘model self-energy’ 

E(w). 

Bottom: Graphical solution of equation 1(22} for the quasipar¬ 
ticle energies uik with the ‘model self-energy’ E(w). 



perature the peak G which opens the insulating gap de¬ 
creases and shifts to the lower or upper edge of the gap 
whereas the value of ry-j- — 774 , increases. NiO thus changes 
smoothly from a Mott-insulator to a spin-density-wave 
insulator. 

Next, we discuss the origin of the dispersionless bands 
observed in the ARPES spectra of both Shen et al.[ 81] 
and Kuhlenbeck et al.\ 82]]. For the sake of illustration 
we consider a ‘model self-energy’ obtained by arbitrar¬ 
ily choosing a few cr,; and Q in (I2T1) . Figure fl3l shows the 
imaginary part of the resulting £(w) for u slightly off the 
real axis (top) and £(w) for real oj (bottom). The bottom 
part again shows the density of states of the noninteract¬ 
ing band e k as well as the lines to — e_ and to — e+. For 
real w £(w) takes any value in [— 00 , 00 ] precisely once 
in any interval [£j, Ci+i] so that the line ui — e k intersects 
£(w) once for each e k - The point of intersection thereby 
is between between those of the lines uj — e_ and to — e+. 
This shows that in between any two successive poles of 
the self-energy there is one complete quasiparticle band. 
If a given pole has a small cr,;, however, £(w) drops almost 


vertically near the corresponding Q , so that the width of 
the respective band becomes small. Replacing 


E 


Qj 

u - Ci 


~^C + 


go 

w - Co 


in the neighborhood of such a pole Co, the resulting dis¬ 
persion and quasiparticle weight Z = (1 — §§) are 


w/c 


Zk 


c , a ° 

Co+ Co-C-,*' 
(Co - c - e k y ■ 


Therefore, unless the denominator happens to cross zero 
near Co this results in a band with little dispersion and 
low spectral weight close to Co- Whether the band is on 
the high or low energy side of Co depends on the sign of 
the denominator Co — C — e k - Figure [TT] compares partial 
ARPES spectra along (1,0,0), where the sum in (l20l) is 
restricted to either e s -like or f 2 g-like Ni 3d orbitals and 
the respective self-energies. Although the situation in 
NiO is more complicated due to the multi-band situation 
and the hybridization with the O 2 p bands it is quite 
obvious how the various dispersionless bands can be as¬ 
sociated with poles of the self-energy. In the case of NiO 
these poles describe the multiplet splitting of the final 
state of the photoemission process, i.e. mainly the Ni 3+ 
ion. In the absence of the Coulomb interaction a single 
Ni 3d shell would have eigenstates obtained by distribut¬ 
ing the electrons over the e g and t 2 g levels and the single 
particle spectral function A(u>) would have few peaks cor¬ 
responding to the energies of these CEF levels. The con¬ 
siderably larger number of CEF-split multiplet states in 
the presence of Coulomb interaction - as given e.g. in the 
Tanabe-Sugano diagrams - increases the number of peaks 
in A(ui) and the interacting peak structure is generated 
by the poles of the self-energy of the Ni 3d electrons in 
exactly the same way as in Figure (T3l In the solid these 
poles of the self-energy then generate the dispersionless 
bands observed in ARPES as discussed above. In that 
sense one can literally see the dispersionless self-energy 
of the Ni 3d electrons directly in the experimental data 
of Shen et al. and Kuhlenbeck et al.. 


VII. CONCLUSION 

In summary, the Variational Cluster Approximation 
proposed by Potthoff allows to combine the classic field 
theroretical work of Luttinger and Ward with the very 
successful cluster method due to Fujimori and Minami re¬ 
sulting in an efficient band structure method for strongly 
correlated electron systems. Since the VC A is based on 
exact diagonalization which is free from the minus-sign 
problem it allows to take into account the full Coulomb 
interaction in the TM 3d-shell which is known to be 
crucial for reproducing the correct multiplet structure 
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angle-integrated density of states is very similar for the 
paramagnetic and antiferromagnetic phase. The angle- 
integrated spectrum and band structure in the antifer¬ 
romagnetic phase agree well with experiment, whereby 
the band structure shows a considerable number of both 
dispersive and dispersionless bands and again shows the 
massive impact of the strong correlations in NiO in that it 
differs strongly from the band structure obtained within 
DFT. 

Acknowledgement: The author would like to thank F. 
Hardy and K. Grube for help with the specific heat data 
and R. Heid for providing the k-meshs for Brillouin zone 
integration. 


FIG. 14: Comparison of partial ARPES spectra and respec¬ 
tive self-energy for e B orbitals (top) and t 2 g orbitals (bottom). 
The two different self-energies refer to the two spin directions. 


and for obtaining agreement with experiment for angle- 
integrated valence band photoemission [2p - [29f and X-ray 
absorption j30i - |40j . As might have been expected on the 
basis of the success of the cluster method in describing 
these spectroscopies, the multiplet structure turns out 
to be important also for reproducing the experimental 
valence band structure as observed in ARPES in that 
it produces a number of nearly dispersionless bands ob¬ 
served there. The VCA moreover delivers an estimate for 
the Grand Potential and - as demonstrated above - allows 
the discussion of thermodynamics and phase transitions. 
It thereby gives a unified description for a wide variety of 
experimental quantities which probe energy scales from 
the meV range up to ~ 10 eV. 

In the case of NiO, using realistic values of the Hubbard- 
17 and charge transfer energy A - as demonstrated by 
the position of the satellite and the magnitude of the 
insulating gap - and a moderately adjusted value of the 
Slater-Koster parameter (pda) (increased by 10% as com¬ 
pared to the LDA band structure estimate) the theoret¬ 
ical Neel temperature is 481 Kelvin (experimental value: 
523 Kelvin). The behaviour near T/v is consistent with 
a 2 nd order phase transition in a local-moment system, 
with quite accurate Landau behaviour of the free energy 
and ordered moment below T/v and a Curie-Weiss sus¬ 
ceptibility above T(v- Consistent with experiment the 


VIII. APPENDIX A 

In this Appendix we discuss an unphysical solution 
which appears when 7 parameters are varied. Thereby 
all 4 spin-even paramaters were varied and in addition 
to the spin-odd parameters V-(e g ) and e e ,-(t 2 g) which 
are used in the AF-I solution, also the spin-odd part of 
the e g level energy, e e> _(e g ). The upper part of Figure 
[15] shows the temperature dependence of the spin-odd d- 
Level-to-Ligand hopping V-(e g ). We could have chosen 
any other spin-odd parameter but this one is sufficient to 
discuss what is happening. The lower part of the Figure 
shows the fl'(T) = f2(T) — f(T) where f(T) is a second 
order polynomial which has no physical significance and 
was subtracted from fl(T) for the sole purpose of making 
tiny variations of f2 around this relatively strongly vary¬ 
ing but smooth ‘background’ visible. 

As one would expect, V-(e g ) starts to deviate from zero 
at the Neel temperature T/v = 481 K, see point A in the 
upper part of Figure [T5] At T\ « 434 K a bifurcation 
occurs and two new solutions appear, see point B. At 
T 2 ss 409 K another bifurcation occurs and the solution 
starting out from A disappears (point C). The second so¬ 
lution emerging from the bifurcation C can be followed 
up to temperatures far above T/v. As can be seen in the 
bottom part, the solution along A —> C has a higher Q 
than the one extending from C to high temperatures (the 
part above 430 K is omitted for this solution in the bot¬ 
tom part of Figure fl5l to keep the range of H sufficiently 
small). In fact, H for this solution turns out to be even 
lower than for the paramagnetic solution so that we have 
an antiferromagnetic solution which gives the lowest Q 
up to the highest temperatures studied. Moreover, the 
magnetic solution A —» C appearing at T/v - which would 
be consistent with the staggered susceptibility - would be 
an unstable state which never should be realized. One of 
the two solutions emerging from the bifurcation B inter¬ 
sects the unphysical antiferromagnetic solution at ~ 420 
K. 

The bottom part also shows fi'(T) for the solution AF-I 
with 6 varied parameters. Over almost the entire tem¬ 
perature range this solution is very close in energy (the 
deviation is oc 10 -5 eV between 300 K and 420 K and 
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Temperature [K] 



Temperature [K] 


FIG. 15: Top: Spin-odd hopping parameter V-(e g ) for the so¬ 
lution with 7 parameters as a function of temperature. Bot¬ 
tom: Grand Potential as a function of temperature for the 
different solutions with 7 parameters (lines). To make small 
changes visible the function f(T) = — 2.99067e • 10 -7 T 2 + 
1.06805 • 10 ~ a T — 712.08 was subtracted from SI. Also shown 
is Q'(T) for the solution AF-I with 6 parameters (circles). 


even smaller above 440 K) to one of the solutions with 
7 parameters. In fact, it seems to interpolate between 
two branches of solutions with 7 parameters. That the 
‘closeness’ is not restricted to fl can be seen in Table 
El which compares some physical quantities for the two 
different solutions, and from Figure [16] which compares 
the k-integrated spectral densities. Both, the table and 
the Figure, shows an essentially perfect agreement be¬ 
tween the two solutions as far as observable quantities are 
concerned. In contrast to this, the symmetry-breaking 
parameters A^_ which are also listed in the table are 
quite different and in fact substantially larger for the so¬ 
lution with 7 parameters. This shows that the symmetry¬ 
breaking effects of V-(e g ) and ed-(e g ) must cancel to a 
large degree in order to simulate the effect of a signifi¬ 
cantly smaller V- (e g ) alone. Adding ed,-(e g ) to the set of 



AF-I 

7 Pars 

V-(e 9 ) 

0.0338 

0.1021 

e d ,-(e g ) 

0.0000 

0.0414 

ed,-(t 2 g 

0.0171 

0.0198 

SI + fj,N 

-244.043562 

-244.043570 

( H ) 

244.0169 

244.0169 

S/ks 

0.7734 

0.7736 

(Ho) 

-457.6312 

-457.6322 

(Hi) 

213.6143 

213.6153 

(H pd ) 

-3.4749 

-3.4750 

(Hpp) 

-0.1822 

-0.1822 

(H dd ) 

-0.0032 

-0.0032 

K„,t) 

0.5325 

0.5328 

(ne g ,l) 

1.6604 

1.6601 

(nd) 

8.1928 

8.1929 

(m s ) 

-1.1279 

-1.1273 


TABLE II: Comparison of the AF-I solution with 6 parame¬ 
ters and the solution for 7 parameters at 400 K. All energies 
in eV. 


parameters to be varied is superfluous. For completeness 
we note that keeping ed-(t 2g ) = 0 so that V-(e g ) re¬ 
mains as the only spin-odd parameter to be varied, gives 
results which are almost identical to those for the solu¬ 
tion AF-I. From the above it looks very much as if already 
with 6 parameters the solution is converged with respect 
to the number of parameters and that adding an addi¬ 
tional spin-odd parameter results in no more significant 
changes to obervables but creates a new branch of un¬ 
physical solutions. 

We now discuss the artificial antiferromagnetic solution 
at high temperatures. Table IIIII compares some observ¬ 
ables for this solution and the paramagnetic one at 600 K. 
As already mentioned, the unphysical solution has lower 
fl than the paramagnetic phase and in addition also a 
higher entropy. There are no large differences in the 
various ground state expectation values. Although the 
symmetry breaking parameters V-{e g ) and e d ,-(e g ) are 
substantially larger than those for the antiferromagnetic 
solutions at 400 K given in Table [TT] the ordered moment 
is much smaller. This shows again that the symmetry 
breaking effect of the different spin-odd parameters can¬ 
cels almost completely in this solution. 

The above example shows that including too many sym¬ 
metry breaking parameters into the subset of variational 
parameters can lead to unphysical solutions. Here it is 
interesting to note that Kozik et al. recently reported un¬ 
physical solutions obtained by DMFT for the self-energy 


of a simple Hubbard-dimer[85|. A tendency to produce 


additional unphysical solutions may be a general fea¬ 
ture of schemes which aim at computing the self-energy. 
Within the VCA the simplest solution to this problem 
would be to simply restrict the number of parameters to 
a minimum - as was done in the solution AF-II in the 
main text. Clearly this introduces a certain arbitrary- 
ness regarding the subset of parameters to be varied. On 
the other hand the above discussion shows that results of 
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Energy [eV] 


the VCA appear to converge rather well with the number 
of variational parameters. In any way it seems desirable 
to find criteria which allow to identify unphysical solu¬ 
tions or methods to regularize the variational procedure 
so that unphysical solutions are suppressed. 



AF 

Para 

V-{e a ) 

-0.3865 

0.0000 

ed,-{e g ) 

-0.2424 

0.0000 

ed,-(t2g 

-0.0126 

0.0000 

fi + fiN 

244.0621 

-244.0613 

(H) 

-244.0046 

-244.0045 

S/k B 

1.1139 

1.0981 

(Ho) 

-457.3758 

-457.3815 

(Hi) 

213.3712 

213.3770 

( H pd ) 

-3.4373 

-3.4379 

(Hpp) 

-0.1756 

-0.1757 

(H dd ) 

-0.0034 

-0.0034 


1.1138 

1.0945 

( n e g ,i) 

1.0757 

1.0945 

(nd) 

8.1889 

8.1890 

(m a ) 

0.0376 

0.0000 


TABLE III: Comparison of the unphysical AF solution and 
the paramagnetic solution at 600 K. All energies in eV. 


FIG. 16: Comparison of the k-integrated spectral densities for 
the AF-I solution with 6 varied parameters and the solution 
with 7 varied parameters at 400 K. 
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